#The following figures were generated in ArcMap for the article. The following commands enable generating similar figures in R

library(maps)
library(mapdata)
library(maptools)
library(sp)
library(rgdal)

###############
###FIGURE 2###
###############
intersection <- readOGR("/Users/jackpaine/Box Sync/Settlers_violence/Shapefiles/Final shapefiles/","intersection")
mediterranean<-readOGR("/Users/jackpaine/Box Sync/Settlers_violence/Shapefiles/Final shapefiles/","mediterranean")
map("worldHires", xlim=c(-25,55),ylim=c(-35,38))
plot(mediterranean,add=TRUE,col="gray40")
plot(intersection,add=TRUE,col="gray40")

################
###FIGURE D1###
################
mediterranean<-readOGR("/Users/jackpaine/Box Sync/Settlers_violence/Shapefiles/Final shapefiles/","mediterranean")
map("worldHires", xlim=c(-25,55),ylim=c(-35,38))
plot(mediterranean,add=TRUE,col="gray40")

################
###FIGURE D2###
################
rainfall<-readOGR("/Users/jackpaine/Box Sync/Settlers_violence/Shapefiles/Final shapefiles/","high_rainfall_africa")
map("worldHires", xlim=c(-25,55),ylim=c(-35,38))
plot(rainfall,add=TRUE,col="gray40")

################
###FIGURE D3###
################
elevation<-readOGR("/Users/jackpaine/Box Sync/Settlers_violence/Shapefiles/Final shapefiles/","high_elevation_africa")
map("worldHires", xlim=c(-25,55),ylim=c(-35,38))
plot(elevation,add=TRUE,col="gray40")

################
###FIGURE D4###
################
tsi<-readOGR("/Users/jackpaine/Box Sync/Settlers_violence/Shapefiles/Final shapefiles/","low_tsi")
map("worldHires", xlim=c(-25,55),ylim=c(-35,38))
plot(tsi,add=TRUE,col="gray40")